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ABSTRACT 

The Galaxy's stellar halo seems to be a tangle of disrupted systems that have been 
tidally stretched out into streams. Each stream approximately delineates an orbit in 
the Galactic force-field. In the first paper in this series we showed that all six phase- 
space coordinates of each point on an orbit can be reconstructed from the orbit's path 
across the sky and measurements of the line-of-sight velocity along the orbit. In this 
paper we complement this finding by showing that the orbit can also be reconstructed 
if we know proper motions along the orbit rather than the radial velocities. We also 
show that accurate proper motions of stream stars would enable distances to be de- 
termined to points on the stream that are independent of any assumption about the 
Galaxy's gravitational potential. Such "Galactic parallaxes" would be as fundamental 
as conventional trigonometric parallaxes, but measureable to distances ^ 70 times 
further. 
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1 INTRODUCTION 

The precision photometry for millions of faint stars ob- 
served by the Sloan Digital Sky Survey (SDSS) led to 
the discovery of numerous tidal streams of halo stars. 
Many of these streams are certainly of tidal origin be- 
cause t he progenitor h a s bee n seen ijOdenkirchen et al] 



20021: iMaiewski et al] 12004 iBelokuroy et al.| |2006~ 



Fellhauer et al.l l2007l : iGrillmair fc JohnsonI l2006l l but 



some other cases the progenitor i s unknown and may no 
longer be extant ("Grillmaij l2006l : iBelokurov et all l2006bl : 
IGrillmair fc Dionatos 2006i) . It has long been recognised 
that streams provide an important diagnostic of the still 
uncertain Galactic gravitational field by virtue of the 
closeness with which a thin stream approximates an orbit 
l|johnston et all liggd ). However, the traditional way of 
exploiting this connection, which is to search for orbits that 
are consistent with the data, has yielded fewer convincing 
fits to the data than one might have expected, and in any 
given case it is not clear why a better-fitting orbit has not 
been found. 

Recently it was realised that given a Galactic poten- 
tial $ and line-of-sight velocities along a stream, one can 
uniquely solve for the six phase-space coordinates that 
points on the stream mus t have if they are to trace an or- 
bit in the given potential (|jin fc Lvnden-Bellll2007l : iBinnevI 
I2OO8I . hereafter Paper I). If the wrong gravitational potential 
is used in the reconstruction, the recovered phase-space co- 
ordinates will in general be inconsistent with conservation of 



energy (Paper I) and will violate the tangential component 
of the equations of motion (Eyre & Binney 2009, hereafter 
Paper II). Hence the reconstruction technique provides a 
powerful diagnostic of the gravitational potential, and once 
the potential has been determined, it will provide distances 
to stars that He on streams that are as absolute as trigono- 
metric parallaxes but very much more precise than will be 
possible for such distant objects in the foreseeable future 
(Paper I). 

The main obstacles to attainment of these exciting goals 
are (a) the fact that streams differ slightly but significantly 
from orbits, and (b) a lack of reliable line-of-sight velocities 
along streams. Paper II addresses problem (a). This paper 
addresses problem (b) by showing that proper motions may 
be employed instead of line-of-sight velocities. 

Many of the most promising streams have distances in 
the range 10 — 50 kpc, so their distance moduli are 15 — 18.5 
and their solar-type stars have apparent magnitudes in the 
range J ~ 19 — 22.5. Perhaps the clos est streams of interest 
are the GD-1 and An ticentre streams (|Grillmair fc Dionatod 
l2006l : lGrillmai3l2006l ). which are onl y 10 kpc distant. Con- 
sequently at r < 19 IKoposov etaP (|2009l ) were able to ob- 
tain velociti es for 24 turnoff stars in the GD-1 stream, while 
at g < 20 IGrillmair et al.l (|2008h measured velocities for 
~ 20 stream stars in each of two fields. The situation re- 
garding velocities of stars i n the more distan t Orph an stream 
is much less satisfactory - IBelokurov et al] ((20o3) conclude 
that indications of the line-of-sight velocity of the Orphan 
stream "are suggestive rather than conclusive". Even with 
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an 8 m telescope it is extremely challenging to measure the 
line-of-sight velocities to a few kms^^ of significant num- 
bers of main-sequence stars at distances in excess of 20 kpc. 
Consequently, the strategy generally adopted with more dis- 
tant streams is to identify giant stars that probably belong 
to the stream and meas ure their velocities. For example, 
lOdenkirchen et alj (|2009l ) used the VLT to measure the line- 
of-sight velocities of 74 giant stars with i < 18.4 in the 
region of the Pal 5 stream and concluded that only 17 of 
these stars were stream members; because the stream is de- 
fined by main-sequence stars, not giants, one cannot be sure 
that a giant observed in the direction of the stream is not a 
foreground or background object. Moreover, the number of 
main-sequence stars in a length of stream is large compared 
to the number of giants, so there is much greater scope for 
beating down random errors if main-sequence stars can be 
used. 

Since streams are identified from the photometry of in- 
dividual main-sequence stars, it is in principle possible to 
measure proper motions for all the stars that define the 
strea m. Such work is a lready possible with the SDSS sur- 
vey (|Munn et al.l |2004| ) , and work of significantly greater 
precision will be possible with the Pan-Starrs survey, which 
is currently getting underway. In this paper we show that 
orbit reconstruction is possible given proper motions rather 
than line-of-sight velocities. 

In Section 2 we describe the algorithm. Section 3 in- 
troduces the concept of Galactic parallaxes which arises in 
connection with the algorithm. Section 4 reports tests of the 
algorithm. Section 5 sums up and looks ahead. 



2 THE ALGORITHM 

We work in the inertial coordinate system in which the 
Galactic centre is at rest. Let ro be the position vector of the 
Sun, r that of a star in the stream and let s be the vector 
from the Sun to the star: 



ds. 



ro. 



(1) 



where s is the direction from the Sun to the star. We must 
distinguish between the derivative Ds/Df that takes into 
account the velocity vo of the Sun, and the derivative ds/dt 
that does not; the latter is tangent to the stream and is 
the proper motion that would be measured by an observer 
who is stationary at the current location of the Sun, while 
the former is the observable proper motion of a star and 
has a component perpendicular to the stream. If u is the 
angle along the stream from some fiducial point, then we 
may write 

Ds J. ds . „ 

; — = "P' (2) 



Dt 



dt 



where ^ and t are the magnitude and direction of the mea- 
sured proper motion, while ii and p are the magnitude and 
direction of the motion along the stream. 

The space velocities measured by an observer moving 
with the Sun and one stationary at the Sun's location are 
related by 



Ds ds 

m ^ d^"^°- 

With equations ((^J we can therefore write 



(3) 



-— s -|- sat = — s -I- sup - 
Dt dt ^ 



Vo, 



(4) 



where Ds/Dt is the spectroscopically measured heliocentric 
velocity and ds/dt = is the projection along the line 
of sight of the star's velocity with respect to the Galactic 
centre. Equating components in the plane of the sky, we 
have 



s/it — sup — (vo — s ■ Vo s) = sup — Vo±, 



(5) 



where vo_l is the component of the Sun's velocity perpen- 
dicular to the line of sight. This equation has just two un- 
knowns, ii and s, and we can in principle solve for both 
through 



up = fit + 



vo± 



(6) 



Specifically, since both t and p can in principle be deduced 
from the observations and vo± may be presumed known, we 
could determine s such that the right side is parallel to p, 
and then read off u from the magnitude of the right side. 
Unfortunately, the uncertainty in the direction t is likely to 
be significant. We therefore eliminate it by squaring up. 



(su) - 2vo± ■ psit-h |vo±| - (s/i) =0. 

The roots of this quadratic equation for sii are 



sii = vo± ■ P ± V (vo± ■ p)^ - |vo±| + (s/i)2 



(7) 



(8) 



Using this expression to eliminate ii from equation ®, we 
find 



(voj 



|vox|^ + (s/i)2 = /it ■ p. 



(9) 



Since jj, is inherently positive, the sign in front of the radical 
in equation ((Sjl must be chosen to agree with the sign of 
t-p. Even though the directions of individual proper motions 
may be uncertain, it should be possible to decide whether 
they are on average opposed to the direction of travel along 
the stream. Equation ([8]l makes u into a function of s and 
quantities that can be determined from the observations. 

Now let F(r) be the Galaxy's gravitational acceleration 
(F = — V"!"). We recall that when we resolve the star's equa- 
tion of motion along the line of sight to the Sun, we obtain 
(Paper I) 



d?;|| 
'df 



du| 
du 



" s 



(10) 



where the subscripts || and _L denote components along and 
perpendicular to s, respectively. 

Since ii = v±/s and equation © makes it a known 
function of s, with equation (I10|l we can now write down 
a system of three differential equations for the unknowns 
along the stream: 

ds _ 

du u 
d«|| F|[ + sii^ 
du ii 

- i 

du u 

To integrate these equations along the stream we need initial 
conditions for s, uy and t. We can trivially set f = at a 
fiducial point of the stream and guess a value so for s at 
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that point. Then we can compute the initial value of iiy as 
follows. 

We write 
_ dv , d , ^ . , , . 



ds , 
" '^ii ■ p + v\\u + su 
" at " 



(12) 



— 2ii||U + sil, 

where we have used equation @ to eliminate ds/dt. Thus 



d(su) . d(su) 
= u- 



dt 



du 



(13) 



The left side of this equation is obtained by explicitly dif- 
ferentiating equation 



d(su) 1 / 

where we have defined, 

dvo± , , dp 

a = — P + vox • — , 

dw du 



dvf) 



du 



2dfl , S^i^Uy 

- MS — H — 

au u 



.(14) 



(vo± • p)^ - |voj 
7 = Vox ■ p. 



(15) 



Equation p4p is linear in and with equation (|13p it readily 
yields 



■ a/3 - a7 + vqx 



d/i 2 

~A 

au au 
I {l3 + sii^/u) 



(16) 



Once the distance so to the fiducial point has been chosen, 
the right side of this equation can be evaluated from the 
data because differentiating vox = vo — s ■ vo s along the 
stream yields 

ds 



dvoi 



ds 



Vq s - 



s ■ vo — 

au 

S ■ Vo p. 



du du 

= -p ■ Vq S 

Hence 
dvox . 

— p = -s ■ vo, 

au 

and 

dvox ,^ s 

Vox ■ — ; = -(s ■ Vo) voj 

du 



(17) 
(18) 

(19) 



Thus the initial conditions required for the integration of 
equations (lllf) follow once so has been chosen. The solution 
to these equations completes the information required to 
assign a full six-dimensional phase-space position to every 
point on the stream. 



3 GALACTIC PARALLAX 

If we knew vq accurately, we could determine the distance to 
a stationary star from the magnitude of its proper motion, 
which would be entirely due to the Sun's motion. Similarly, 
in so far as we can argue that a star in a stream has no 
velocity perpendicular to the stream, we can determine the 
distance to the star by attributing to the Sun's motion the 
star's proper motion perpendicular to the stream - equation 
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Figure 1. Right panel: the projection onto the sky of a nu- 
merically integrated orbit in a Myamoto-Nagai potential with 
b/a = 0.2. Left panel: the proper motion as a function of b. 
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Figure 2. The log to base 10 of the rms variation in the energy 
when the orbit shown in Fig. ^ is reconstructed from an assumed 
fiducial distance s rather than its true value sq- The triangles 
show results obtained when the reconstruction employs the true 
potential, while the squares and pentagons are for less-flattened 
potential. 



embodies this idea mathematically. Distances obtained 
in this way without the use of dynamics would have the same 
logical status as conventional trigonometric parallaxes, and 
might be called "Galactic parallaxes" . For given astrometric 
precision Galactic parallaxes could be accurately measured 
to much greater distances than conventional parallaxes be- 
cause in three years the Sun moves ~ 140 AU around the 
Galaxy, leading to a change in the position of an object 
that is ~ 70 times larger than the corresponding conven- 
tional parallax angle. Consequently equipment such as Gaia 
that can measure the conventional parallaxes of sources at 
distances of order 10 kpc could measure Galactic parallaxes 
out to 700 kpc, i.e., as far as the Andromeda galaxy. Un- 
fortunately, before this method could be applied to streams 
within M31, one would have to determine the velocity of 
the Sun relative to the bary centre of M31, which we are 
not likely to know better than we know the distance to the 
centre of M31. 

By measuring the Galactic parallax to each point along 
a stream, the stream's three-dimensional trajectory could be 
determined without any knowledge of the Galaxy's gravita- 
tional potential. The requirement that this trajectory be an 
orbit must constrain the potential rather tightly. 

In the previous section we have chosen to sacrifice some 
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of the diagnostic power of equation ® by eliminating t on 
the grounds that it will be hard to measure. With t elim- 
inated, the distance can only be recovered by adopting a 
trial gravitational potential and searching over the fiducial 
distance sq. 

For sufficiently small s, the radical in equation ((Sj be- 
comes imaginary, so there is a lower bound on the values 
of So that should be considered: below this bound the val- 
ues of s obtained by solving equations will somewhere 
approach the value at which the radical in equation ((Hjl be- 
comes imaginary. This event signals that the measured value 
of /i is too small to be consistent with the reffex of the Sun's 
velocity at the proposed distance. Thus the kinematics of the 
problem imposes a lower bound on so- There is no similar 
kinematic upper bound on so because the radical in equation 
© is real for all large s. 



4 TESTS 

We have tested the ability of the algorithm to reconstruct 
orbits in the same Miyamoto-Nagai (1975) potential that 
was used in Paper I, namely 

C M 

HR, = — (20) 

with q = 0.2. The numerical procedures were essentially 
the same as those described in Paper I except for one sig- 
nificant item: the fiducial point was placed in the middle of 
the stream rather than at one of its ends, and equations (|lip 
were integrated in both directions away from this point. This 
modification is advantageous as the data constrain deriva- 
tives of the observables much more tightly at the middle of 
the stream than at its ends. 

The right panel of Fig.[T]shows the track over the sky of 
an orbit that starts at galactic latitude 6 = — 5° and a dis- 
tance So = 15a as viewed from the Sun, which is on a circular 
orbit at Rq = 8a. The left panel shows the magnitude of the 
proper motion along the orbit. The points in Fig.[2]show the 
rms variation in the energy of orbits that are reconstructed 
from the 42 data points shown in Fig.[l]for various assumed 
distances s to the fiducial point and three trial potentials: 
the true potential, which has scale-length ratio q — 0.2 and 
two less fiattened potentials. The rms variation in E has a 
sharp minimum at the true distance when the true poten- 
tial is used, and higher minima when the wrong potential is 
used. The exquisite precision with which the distance can be 
determined from this plot is remarkable: better than three 
parts in 10"*. 

In Fig. [2] the smallest value reached by the rms energy 
variation is significantly higher than the corresponding fig- 
ure for the same orbit when line-of-sight velocities are used 
(Paper I). However, for other orbits smaller values of the rms 
variation in energy are obtained from proper-motion data; 
not surprisingly the data with the greatest diagnostic power 
varies with the nature of the orbit. 



5 CONCLUSIONS 

We have complemented the work of Paper I by showing that 
when proper motions can be measured along a section of a 



single orbit, the full phase-space coordinates for the orbit 
can be reconstructed as readily as is the case when line-of- 
sight velocities have been measured. However, we have also 
recovered a more powerful result: if the direction as well 
as the magnitude of the proper motions can be accurately 
measured, the three-dimensional geometry of the orbit can 
be recovered without assuming anything about the Galaxy's 
gravitational potential. This reconstruction is possible be- 
cause the Sun's velocity must be responsible for the motion 
of stars perpendicular to the orbit, so from proper motion 
perpendicular to an orbit distances can be inferred that are 
as fundamental as conventional trigonometric parallaxes. 

The extension of the work of Paper I to proper motions 
is useful because proper motions can be readily measured for 
the main-sequence stars that define the stream. These stars 
are typically so faint (7 > 20) that it is hard to obtain line- 
of-sight velocities of the requisite accuracy for large numbers 
of them. 

In Paper II we show that even when line-of-sight veloci- 
ties are available, a combination of observational errors and 
the fact that a tidal stream does not strictly follow an orbit 
make the task of determining whether a given potential is 
compatible with an observed stream complex. It proves nec- 
essary to search a multi-dimensional space of phase-space 
tracks that are compatible with the data for ones that are 
dynamically consistent orbits. A comparable search would 
be required when the data included proper motions rather 
than line-of-sight velocities. In a future paper we will extend 
to proper motions the techniques for handling this problem 
that are developed for line-of-sight velocities in Paper II. 

Another interesting avenue is to determine the solar 
motion by treating it as an unknown when reconstructing 
orbits for several streams, using either line-of-sight veloc- 
ities or proper motions, and choosing the value which al- 
lows a consistent int erpretation of these streams. Recently 
IKodosov et al.l (j2009l ) used a combination of line-of-sight ve- 
locities and proper motions for the GD-1 stream to constrain 
vo. Stronger constraints should be attainable by modelling 
several streams simultaneously. 
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